arXiv:1507.05804v2 [math.PR] 4 Aug 2015 


Spatial birth-and-death Markov dynamics of finite particle 

systems 


Viktor Bezborodov* 

Bielefeld University, Faculty of Mathematics 


August 5, 2015 


Abstract 

Spatial birth-and-death processes with time dependent rates are obtained as solutions 
to certain stochastic equations. Existence, uniqueness, uniqueness in law and the Markov 
property of unique solutions are proven when the integral of the birth rate over R d grows 
not faster than linearly with the number of points of the system. Martingale properties of 
the constructed process provide a formal connection to the heuristic generator. 

We also study pathwise behavior of aggregation model. We estimate the probability of 
extinction and the speed of growth of the number of particles of the process conditioned on 
non-extinction. 


1 Introduction 

We consider spatial birth-and-death processes with time dependent rates. At each moment of 
time the system is represented as a finite collection of motionless points in M d . We interpret 
the points as particles, or individuals. Existing particles may die and new particles may appear. 
Each particle is characterized by its location. 

The state space of a spatial birth-and-death Markov process on M d with finite number of 
points is the space of finite subsets of M d , 

r 0 (M d ) = {r, C M d : 17/| < oo}, 

where \rj\ is the number of points of y. 
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Denote by ^(M d ) the Borel c-algebra on M. The evolution of the spatial birth-and-death 
process on M d admits the following description. Two measurable functions characterize the 
development in time, the birth rate b : M d x M + x To(M d ) —>• [0, oo) and the death rate d : 
M d x M_|_ x To(M d ) —>■ [0, oo). If the system is in state rj £ To at time t, then the probability that 
a new particle appears (a “birth”) in a bounded set B £ &(M. d ) over time interval [f; t + At] is 


At J b(x, t, rf)dx + o(At), 

B 

the probability that a particle x £ rj is deleted from the configuration (a “death”) over time 
interval [t]t + At] is 

d(x, t, 77 ) At + o(At), 

and no two events happen simultaneously. By an event we mean a birth or a death. Us¬ 
ing a slightly different terminology, we can say that the rate at which a birth occurs in B is 
f B b(x,t,rj)dx, the rate at which a particle x £ i] dies is d(x,t, rj), and no two events happen at 
the same time. 

Such processes, in which the birth and death rates depend on the spatial structure of the 
system as opposed to classical Z + -valued birth-and-death processes (see e.g. [Har63, Page 116], 
[AN72, Page 109], and references therein), were first studied by Preston [Pre75]. A heuristic 
description similar to that above appeared already there. Our description resembles the one in 
[GK06]. 

We say that the rates b and d, or the corresponding birth-and-death process, are time- 
homogeneous if b and d do not depend on time. By abuse of notation we write in this case 
b(x,s,rj ) = b(x,rj), d(x,s,ij) = d(x,rj). The (heuristic) generator of a time-homogeneous spatial 
birth-and-death process should be of the form 


LF(rj) = / b(x, ri)[F(r] U x) — F(rj)\dx + d(x , rj){F{j] \ x) — F(rj)), 

xGri 


(1) 


for F in an appropriate domain, where rj U x and rj \ x are shorthands for rj U {x} and rj \ {x}, 
respectively. 

Garcia and Kurtz [GK06] obtained birth-and-death processes as solutions to certain stochas¬ 
tic integral equations. The systems treated there involves an infinite number of particles. In the 
earlier work [LG95] of Garcia another approach was used: birth-and-death processes were ob¬ 
tained as projections of Poisson point processes. A further development of the projection method 
appears in [GK08]. Fournier and Meleard [FM04] used a similar equation for the construction 
of the Bolker-Pacala-Dieckmann Law process with finitely many particles. Following ideas of 


2 


[GK06] and [FM04], we construct the birth-and-death process described above as solution to a 
stochastic equation. 

Holley and Stroock [HS78] constructed the spatial birth-and-death process as a Markov 
family of unique solutions to the corresponding martingale problem. For the most part, they 
consider a process contained in a bounded volume, with bounded birth and death rates. They 
also proved the corresponding result for the nearest neighbor model in M 1 with an infinite number 
of particles. 

Belavkin and Kolokolt.sov [BK03] discuss, among other things, a general structure of a Feller 
semigroup on disjoint unions of Euclidean spaces (see also references therein for the construc¬ 
tion of the Markov processes with a given generator). We note that time-homogeneous birth- 
and-death processes need not have the Feller property. Eibeck and Wagner [EW03] discuss 
convergence of particle systems to limiting kinetic equations. In particular, they construct the 
stochastic process corresponding to the particle system as a minimal jump process, or pure jump 
type Markov process in the terminology of Kallenberg [Kal02], The jump kernel is assumed to 
be locally bounded. 

Kondratiev and Skorokhod [KS06] constructed a contact process in continuum with an in¬ 
finite number of particles. The contact process in continuum can be described as the spatial 
birth-and-death process with 

b(x,rj) = A^a(x - y), d{x,rj) = 1, 

y£il 

where A > 0 and 0 < a € L 1 (R d ). Under some additional assumptions, they showed existence 
of the process for a broad class of initial conditions. 

The scheme proposed by Etheridge and Kurtz [EK14] covers a wide range of interactions 
and applies to discrete and continuous models. Their approach is based on, among other things, 
assigning a certain mark (’level’) to each particle and letting this mark evolve according to 
some law. A critical event, such as birth or death, occurs when the level hit some threshold. 
Shcherbakov and Volkov [SV15] consider the long term behavior of birth-and-death processes 
on a finite graph with constant death rate and the birth rate of a special form. Bezborodov 
[Bezl5] obtains the spatial time-homogeneous birth-and-death process as a unique solution to 
the equation slightly different from the one we use here. Various questions not treated in this 
paper are considered there, for example the possibility of an explosion, continuous dependence 
on initial conditions and related semigroup of operators. 

In the aforementioned references as well as in the present work the system is represented by 
a Markov process. An alternative approach consists in using the concept of statistical dynamics 
that substitutes the notion of a Markov stochastic process. This approach is based on considering 
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evolutions of measures and their correlation functions. For details see e.g. [FKK12], [FKKZ14], 
and references therein. 

Finkelshtein et al. [FKKZ14] consider different aspects of statistical dynamics for the aggre¬ 
gation model. In this model the death rate is given by 

d(x, rj) = exp I - ^ (fix - y) I , 

\ y£v\ x ) 

where <fi is a positive measurable function. For more details see [FKKZ14] and references therein. 
In this paper we present a detailed long time microscopic description of the process. Also, we 
estimate the probability of extinction and the speed of growth of the average number of points. 

The paper is organized as follows. Notation, definitions and results are given in Section 2. 
Proofs are given in Sections 3 and 4, with two auxiliary results located to Section 5. 

2 The set-up and main results 

2.1 Construction and basic properties 

Let A be a locally compact separable metric space (typically A will be a subset of M d ). Denote 
by r(A) the space of locally finite subsets of X 

T(X) = {7 C A | | 7 n K | < oo}, 

also called the space of configurations over X. The space T(A) can be endowed with a cr-field 
generated by the projection maps 


T(A) 9 7 |y n B\ € z + 

where B is an arbitrary Borel subset of A. For more details about the notions introduced 
here see e.g. [DVJ08], [Kal02, Chapter 12] or [KK02], and references therein. Throughout this 
paper r 2 stands for r(M + X M + ). Let 7 r be the distribution of a Poisson random measure on 
(r 2 ,^(r 2 )), with the intensity measure being the Lebesgue measure on X M + . We denote 
by 3S t (fT 2 ) the smallest sub-er-algebra of ^(r 2 ) such that for every A± € ^([0,f]), A 2 € ^(M+) 
the map 

r 2 3 7 i->- y(Ai x A 2 ) G R+ 
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is ^(r 2 )-measurable. Similarly, define ^>t(r 2 ) as the smallest sub-cr-algebra of such 

that for every A\ G &((t, 00 )), A 2 G ^(R+) the map 

T 2 9 7 7(^1 x ^ 2 ) G M+ 


is ^>^(r 2 )-measurable. 

Let 7/o be a (random) initial configuration, and let 7/0 be the point process on R d xl^ obtained 
by associating to each point in 7/0 an independent Poisson point process on R + x R + , with the 
distribution 7 r. That is, for 7/0 = ; set 

OO 

%= ^(xi, 7i)- 

i=l 

Consider the stochastic equation with Poisson noise 

WtiB) — j d [o,b(x,s,ri 3 -)]( u ) d { 

(0,t]xBx[0,oo)xr2 

/ [0,d(x,r,7 7r _)](^)7(rf ? ’ ) ^) = 0}r/ 0 (dx,d7), 

BxT 2 

where (r]t)t>o is a cadlag To-valued solution process, IV is a Poisson point process on R + X R d x 
R_l_ x r 2 , the mean measure of IV is ds x dx x du x n. We require the processes N and 7/0 to be 
independent of each other. Equation (2) is understood in the sense that the equality holds a.s. 
for every bounded B G ^(R d ) and t > 0. 

Remark. In the first integral on the right-hand side of (2) x is the place and s is the time 
of birth of a new particle. This particle is alive as long as f* I[o,d(x,r,Vr-)\( v ) n f( dr i dv ) = 0, where 
{x, s, u, 7) G N. Thus, 7 is the process ’responsible’ for death. In the death term lies the main 
difference to the equation considered by Garcia and Kurtz [GK06]. Adapted to our notation, 
the equation there is of the form 



/' 


I [0,d{x,r,r lr -)}(. v )'y( dr ’ dv ) = ®}N(ds, dx, du, drj) 


( 2 ) 


rit(B) = 1 d [o,b(x,r] a —)] (w)/{ j d(x,rj r -)dv < r}N(ds, dx, du, dr) 


(0,t] xBx [0,oo) x [0,oo) 


+ / d(x,r] r -)dv < r}rjo(dx,dr), 


(3) 


Bx[0,oo) 


where N is a Poisson point process on R + x R d x R + x R + with mean measure ds x dx xdux e~ r dr, 
and 7/0 is obtained from 7/0 by attaching an independent unit exponential to each point. At first 
glance, (2) is more complicated than (3), since the death mechanism requires a whole Poisson 
random measure on [0; oo) 2 instead of just one exponential random variable. However, it is not 
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clear how to a priori define a filtration {&t}t >0 compatible with N and such that a solution to 
(3), if unique, should possess the Markov property with respect to o- 

Conditions on b, d and t]q. Unless stated otherwise, we assume that the birth rate b satisfies 
the following conditions: sublinear growth on the second variable in the sense that 

/ sup b(x, s,rj)dx < c\\v\ + c 2 , (4) 

J s> 0 
R d 

for some constants ci, C 2 > 0, and that b(x, -, rj) and d(x, -,rj) are right-continuous for any x E R d 
and r/ S To. 

We also assume that 


E 177o| < oo. (5) 

We say that N is compatible with a right-continuous complete filtration {^t} if 

N([0,q] x B x C x S) 

is ^-measurable for any q E [0,t],B E 3S(M. d ),C E 38(R+), and 3 E 38 r 2 ), and if 

N((t,t + q'j xB'xC'x 3') 

is independent of for any q' > 0, B' E 38(JH A ), C' E 38 { R+), and 3' E =^>^ 2 ). 

Sometimes we will use the representations 

^ = ^^^(Sq,Xq, u q,7q)’ = ^ , 7<j) ’ 

gGX qej 

where X and are some countable disjoint sets. 

Definition 2.1. A (weak) solution of equation (2) is a triple ((vt)t>o, n), 0 ), 

where 

(i) is a probability space, and 0 is an increasing, right-continuous and 

complete filtration of sub-u-algebras of J^", 

(ii) N is a Poisson point process on R + x R d x R + x T 2 with intensity ds x dx x du x 7 r, 

(iii) ?7o is a random J^o-measurable element in To satisfying (5), 

(iv) the processes N and vo are independent, N is compatible with {^t}t>o, 

(v) (vt)t> 0 is a cadlag To-valued process adapted to 0 , Vt | t _ 0 = Vo, 
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(vi) all integrals in (2) are well-defined, 


E / ds 


J b(x, 


S,Vs-)+ ^2 d ( x i s ^ S - 

X&TJs- 


< oo, t > 0 


and 


(vii) equality (2) holds a.s. for all £ G [0, oo] and all Borel sets B. 
Let 


y? = °{rio,N([0,q\ xBxCxE), 

q G [0 ,t\,B G f(K d ),C£ J(1 + ),S£ 9B t { r 2 )}, 

and let 5^t be the completion of under P. Note that 0 is a right-continuous filtration. 

Definition 2.2. A solution of (2) is called strong if (r]t)t >0 is adapted to (^,£ > 0). 

Remark 2.3. In the definition above we considered solutions as processes indexed by £ € [0, oo). 
The reformulations for the case £ G [0,T] , 0 < T < oo, are straightforward. This remark also 
applies to many of the results below. 

Definition 2.4. We say that pathwise uniqueness holds for equation (2) and an initial distri¬ 
bution v if, whenever the triples ((rjt)t>o,N), (fl,^,P), o) and ({fj t )t>o, N), (Q,^,P), 

({# t } t>0 ) are weak solutions of (2) with P{t]q = fjo} = 1 and Law(r) q) = u, we have P{r]t = 
fjt,t G [ 0 , oo)} = 1 (that is, the processes 77,77 are indistinguishable). 

Definition 2.5. We say that joint uniqueness in law holds for equation (2) with an initial 
distribution v if any two (weak) solutions ((i]t),N) and of (2), Law(rjo ) = Law(r/' 0 ) = n, 

have the same joint distribution: 


Law((rjt), N) = Law{(j]' t ),N'). 

Theorem 2.6. Pathwise uniqueness, strong existence and joint uniqueness in law hold for equa¬ 
tion (2). The unique solution is a Markov process. If b and d are time-homogeneous, then the 
family of push-forward measures {P a ,a G To} defined in Remark 3.3 constitutes a Markov pro¬ 
cess, or a Markov family of probability measures, on Dy 0 [0, 00 ). 

Proof. The statement is a consequence of Proposition 3.2, Remark 3.3 and Proposition 3.8. 
I 11 particular, the Markov property of {P a ,a G To} follows from Corollary 3.7. 

We call the unique solution of (2) (or, sometimes, the corresponding family of measures on 
Dr 0 [0,oo)) a (spatial) birth-and-death Markov process. 
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Remark. For time-homogeneous b and d, the transition probabilities of the embedded 
Markov chain (see e.g. [Kal02, Chapter 12]) of the birth-and-death process are completely 
described by 

Q(i ?, in \ Ml) = ^B + D){ri) ’ xGT? ’ ( 6 ) 

Q(v, Ul u {*}, xeu})= U € ^(R d ), r, € r 0 , 

where (B + £>)(??) = / Rd b(x, g)dx + Y^ x£r , d (x, d)- 

The following two propositions establish a rigorous relation between the unique solution to 
(2) and L defined by (1). To formulate the first of them, let us consider the class ^ of cylindrical 
functions F : To —>• R+ with bounded increments. We say that F has bounded increments if 

sup (F(r/ U {:r}) — F(rj)) < oo. 

??€ ro,a;GK d 

We say that F is cylindrical if for some R = Rp >0 

F(rf) = F(C) whenever r/ n B(od, R) = C H B(oa, R), 

where B(x, R) is the closed ball of radius R around x, and Od is the origin in R d . We recall that 
the filtration > 0} appeared before Definition 2.2. 

Proposition 2.7. Let {j]t)t> o be a weak solution to (2). Then for any F E the process 

t 

F{f}t) ~ J { J b(x,s,r] s _)[F(r] s - U {x}) - F(rj s _)]dx 

o K d (7) 

- d(x,s,rj s -)[F(r] s -\{x}) - F(ri s -)]}ds 

X Cris¬ 
is an > 0} -martingale. In particular, the integral in (7) is well-defined a.s. 

2.2 Aggregation model 

Here we consider a specific time-homogeneous model, which we call an aggregation model. This 
model has a property that the death rate decreases as the number of neighbors grows. We treat 
here the death rate given below in (8), and, in addition to previous assumptions, we require the 
birth rate to grow linearly on the number of points in configuration in the sense (9). We prove 
in Proposition 2.8 that the probability of extinction is small if the initial configuration has many 
points in some fixed Borel set A C R d . Propositions 2.9, 2.10 and Theorem 2.11 describe the 
pathwise behavior of the process. 




Let 


( 8 ) 


d(x, rj) = exp{— ^ ip(x - y)}, 

2/&7 

where p is a nonnegative measurable function. Proposition 3.2 ensures existence and uniqueness 
of solutions, and that the unique solution is a pure jump type Markov process. 

We want to show that, if the initial configuration has m points in some bounded region, then, 
under some assumption on b and the probability of extinction declines at least exponentially 
fast by m. Also, we will give a few statements describing the pace of growth of the number of 
points of the system. 

More specifically, let A be a measurable subset of M d , the birth rate and the initial condition 
r/o satisfy the same condition as in Proposition 3.2, and, besides that, the inequalities 

J b(x, rj)dx > c\rj n A|, ij€T 0 , (9) 

A 

and 

b(x,ri l ) <b(x,i y 2 ), r/ 1 ,!] 2 6 r 0 , rj 1 C r/ 2 , (10) 

hold for some positive c. We assume also that 

inf ip(x - y) > logo, (11) 

x,y£ A 

where a > 1. 

We say that the process (r]t)t> 0 extincts if inf{f > 0 : ry = 0} < oo. This inhmum is called 
the time of extinction. 

Proposition 2.8. Let C > 0. Then there exists mo = mo(C) £ N such that, whenever m > mo, 

P«{(nt ) t>o extincts } < C m 


for all a satisfying \a D A| = m 

Proposition 2.9. For all a € To, 

P a ({\vt n A| ->• oo} U {3t' : Vt > t ', | rj t fl A| = 0}^ = 1. (12) 

Remark. Note that we do not require b(-,0) = 0; if f A b(x, 0)dx > 0, then (12) implies 

Pa{\yt n A| ->• oo} = 1. 

The next proposition is a consequence of the exponentially fast decay of the death rate. 
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Proposition 2.10. With probability 1 only a finite number of deaths inside A occur: 


P a < fit H A| — fit- n A| = —1 for infinitely many different t > 0 > = 0, a: G Tq. 


Theorem 2.11. Let a G IV For P a -almost all uj G F := { lim fit D A| = oo} we have 

£—>•00 

lim inf |,ftn . A| > 0 . 


£—> oo 


0 c£ 


(13) 


Corollary 2.12. For all configurations a with aflA/ 0 , 


inf MhllL > 0 , 


0 c£ 


£>o e 

Remark. If A has a finite volume and the birth rate is given constant within A, that is 


(14) 


b(x, rj) = C 3 > 0, x G A, 


then from the proofs we can conclude that Theorem 2.11 and Corollary 2.12 still hold provided 
that we replace (13) by 

lim inf ^ > 0. (15) 

oo c 3 Vol d (A) 


and (14) by 


inf n A ' > 0. 


(16) 


t>0 C3Vold(A) 

Here C 3 Vol d (A) is the d-dimensional volume of A. These two growth estimates stand in contrast 
to the mesoscopic behavior of the system [FKKZ14]. Theorem 5.3 in [FKKZ14] says that for 
some values of parameters the solution to the mesoscopic equation started from sufficiently 
small initial condition stays bounded. On the contrary, the microscopic system grows whenever 
it survives, and the density always grows. 


3 Proof of Theorem 2.6 and Proposition 2.7 

Let us start with the equation 

Vt( B ) = j 7 [ 0 , 6 (x,s, Vs -)]( U ) N ( ds ’ dx ’ du ’ d 'y) +Vo( B ), (17) 

( 0 ,(]xBx[ 0 ,oo)xr 2 

where b(x,if) := sup b(x,s,£). Note that b satisfies sublinear growth condition (4) if b does. 
s>0,£o? 

This equation is of the type (2), with b being the birth rate, and the zero function being 
the death rate, and all definitions of existence and uniqueness of solutions are applicable here. 
Later a unique solution of (17) will be used as a majorant of a solution to (2). 
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Proposition 3.1. Under assumptions (4) and (5), strong existence and pathwise uniqueness 
hold for equation (17). In particular, the unique solution (fjt.)t>o satisfies 


E\r]t\ < oo, t > 0. 


(18) 


Proof. For to € {/ b(x, go)dx = 0}, set Q = go, a n = oc, n G N. 

R d 

For to G F := { f b(x, go)dx > 0}, we define the sequence of random pairs {(a n , C^)}, where 
R d 


(T n +\ = inf{t > 0 


/ 


f [o,6(x,c CTn )] ( u ) N (ds, dx, du , d'y) > 0} + a n , a 0 = 0, 


{cTn ,CTn +*] X B X [0,do) X T 2 


and 


Co — <r l 0 , Ccr n+ i — C (T n U { z n+l} 


for z n+ 1 = {x G M d : A7({cr^^i} x {x} X [0, b(x, Cct„)] x ^) > 0}. The points z n are uniquely 
determined almost surely on F. Furthermore, <x n+ i > a n a.s., and a n are finite a.s. on F (in 
particular because &(x,Cr„) > b(x,rj 0 )). For to G F, we define Ct = (<r n for t € [a n , cr n +i)- Then 
by induction on n it follows that a n is a stopping time for each n G N, and <f an is D Im¬ 
measurable. By direct substitution we see that (Ct)t>o is a strong solution to (17) on the time 
interval t G [0, lim a n ). Although we have not defined what is a solution, or a strong solution, 

n—>00 

on a random time interval, we do not discuss it here. Instead we are going to show that 


lim a n = 00 a.s. (19) 

n—> 00 

This relation is evidently true on the complement of F. If P{F) = 0, then (19) is proven. 

If P(F) > 0, define a probability measure on F, Q(A) = 4 G / :=#TlF, and define 

= F. 

The process N is independent of F, therefore it is a Poisson point process on the probability 
space (F, J?, Q) with the same intensity, compatible with {^t}t> o- From now on and until other 
is specified, we work on the filtered probability space (F, J?,{j?t}t>o,Q)- We use the same 
symbols for random processes and random variables, having in mind that we consider their 
restrictions to F. 

The process (Ct)te[o, lim a n ) has the Markov property, because the process N has the strong 

n—t 00 

Markov property and independent increments, see Appendix. Indeed, conditioning on J’ an , 

fb(x,(a n )dx 

E [h^ n+ i=Ca n Ux for some xeB} I ^<J n ] = ^ ’ 

Rd 
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thus the chain { (a n }nez + is a Markov chain, and, given {C <r„}nez + , &n+i ~ &n are distributed 
exponentially: 

E{I{a n+1 -a n >a} I {Cx„}neZ+} = exp{-a / b{x,(ajdx}. 

R d 

Therefore, the random variables = ( a n — a n -i)(f b(x. C, an )dx ) constitute a sequence of inde- 

R d 

pendent random variables exponentially distributed with parameter 1, independent of {Cr„}ng.z + - 
Theorem 12.18 in [Kal02] implies that (Ct)tg[o, lim a„) is a pure jump type Markov process. 

n—too 

The jump rate of 0, lim a n ) is given by 

n—> oo 

c(a ) = / b(x, a)dx. 

R d 


Condition (4) implies that c(a) < ci|a| +C 2 - Consequently, 


c((a n ) < CllCrJ +C 2 = Cl |Co| +Cin + C 2 . 


We see that e ^ 1 j = oo a.s., hence Proposition 12.19 in [Kal02] implies that —)• oo. 

Now we return again to our initial probability space (0, {^}t>o, P)- We have proved 

the existence of a strong solution. The uniqueness follows by induction on jumps of the process. 
Namely, let (Ct)t>o be another solution of (17). From (vii) of Definition 2.1 and the equality 


/ 

(0,cri)xM d x [0,oo] XT 2 

it follows that P{( has a birth before 

/ 

{cti } xR d x [0,oo] xT 2 


J [o,b(x, Vo )] ( u)N(ds , dx, du , d'y) = 0, 
cti} = 0. At the same time, the equality 
I [o.b(x,r/o)] { u ) N {d s , dx, du, d'y) = 1, 


which holds a.s., yields that C has a birth at the moment a i, and in the same point of space at 
that. Therefore, C coincides with C up to o\ a.s. Similar reasoning shows that they coincide up 
to a n a.s., and, since a n —y oo a.s., 


P\Ct = Ct for all t > 0} = 1. 

Thus, pathwise uniqueness holds. 

Now we turn our attention to (18). We can write 
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( 20 ) 


Ktl = M + 7 {l&l - M ^ n ) 

n= 1 
oo 

= \m\ + '^2 I { (J n < t}. 

n =1 


Since <J n = J2 


7 i 


,_i f KxXvJdx 

Rd 


we have 


fan < i} = {X] 


7i 


^ / b(x,(*i)dx 


< 


*>c{x; 


7i 


Ci|C CTi | +C 2 


<*} 


ciE 


7i 


— < f} = {Z t - Z 0 > n}, 


(ci + c 2 )(|t? 0 | +i) 

where (Zt) is the Yule process, i.e. the birth process on Z + with transition rates Pk,k+ 1 = 
(ci + c 2 )fc, pkj = 0, l k + 1, see, e.g., [AN72, Chapter 3, Section 5]. Here (Z t ) is defined as 
follows: Zt — Zq = n when 


E 


7 i 


71+1 

<‘<E 


7* 


—( (ci + c 2 )(|?7o| + i) “ (ci + c 2 )(|r/o| + i) ’ 


and Zo = \rjo\. Thus, we have \ft\ < Zt a.s., hence E\ft\ < EZf < oo. The constructed solution 
is strong. □ 

Proposition 3.2. Under assumptions (4)-(5), pathwise uniqueness and strong existence hold 
for equation (2). The unique solution (rjt) satisfies 


E\r)t | < oo, t > 0 . 


( 21 ) 


Proof. Let us define stopping times with respect to {+, t > 0 }, 0 = 9q < 9i < 0 2 < @3 < ..., 
and the sequence of (random) configurations fae+eN as follows: as long as 

9 n +l = #n+l ^ ^n+l + < CO; 


where 


9 b n+ i = inf{7 > 0 : J +, b{x,s, Ve J) {u) N (ds, dx, du, dr/) > 0 }, 

(On ,0n+t] x M d x [0,oo) X F 2 


K+l = inf{7 > 0 : ^ / +,d(z i ,r V?9n )+b q(dr, dv) > 0}, 

q^XU^J d (0n,0n~\-t] X [0,oo) 

xq&lefi 
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we set rie n+l = Ve n U {z n+ 1 } if O b n+ 1 < 9 * +1 , where {z n+1 } = {z € M d : N({9 n + 6 b n+1 } x 
{z} x M+ x r 2 ) > 0}; r]g n+1 = r]g n \ {z n+1 } if 9 b n+1 > fl£ +1 , where {z n+1 } = {x q € ve n ■ 
7 q ({0 n + 0n+i\ x > 0 }; the configuration rjg 0 = r]o is the initial condition of ( 2 ), rjt = r]g n 
for t £ [6 n ,6 n+ 1). Note that 

p {°n +1 = #n+1 I min{^ +1 ,^ +1 } < 00} = 0 , 

the points z n are a.s. uniquely determined, and 

p {zn+ 1 G r)g n | d b n+l < 6^ +1 } = 0. 


If for some n 


0 n -\-i — 00, 


we set 6 n+k = 00, k £ N, and Vt = Ve„, t > 9 n . 

Random variables 9 n ,n £ N, are stopping times with respect to the filtration > 0 }. 

Using the strong Markov property of a Poisson point process, we see that, on { 9 n < 00}, the 
conditional distribution of given ^g n is 


rdn+p 

P{Q b n+ 1 >p} = exp{- / b(x, s, r}g n )ds}, 
J Q n 

and the conditional distribution of 9 ^, 1 given is 


rdn+p 

P{ 9 n +1 > P} = exp{- / d(x, a, 
J 0 n 


In particular, 9 b n , 9 ^ > 0 , n € N. 
We are going to show that a.s. 


9 n —> 00, n —>• 00. (22) 

Denote by the moment of the k- th birth. It is sufficient to show that 9 ' k —>• 00, k —>• 00, 
because only finitely many deaths may occur between any two births, since there are only finitely 
particles. By induction on k' we can see that C {erjjigN, where 07 are the moments of 

births of (Vt)t> 0, the solution of ( 17 ), and rjt C rj t for all t £ [ 0 , lim n 0 n ). For instance, let us 
show that (rj t )t >0 has a birth at 9 [. We have rjgi _ Dfj 0 = 770, and r/g^ _ C Vt |t=o= Vo, hence for 
all x £ M d 


&(®>» 7 fl'-) > b(x, ve[-) > Kx,ve[-) 
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The latter implies that at time moment 8[ a birth occurs for the process (rj t )t >o in the same 
point. Hence, rjg / C %/, and we can go on. Since cifc —>• oo as k —> oo, we also have 6' k —> oo, 
and therefore 6 n oo, n —> oo. 

Let us now prove the inequality from item (vi) of Definition 2.1, 


Z 

E J ds j b(x, s, r/ s -) + ^ d(x,s,%_) 


0 Rd 


x&Vs- 


< oo, t > 0. 


(23) 


Denote the number of births and deaths before f by and dt respectively, i.e. 


bt = #{s ■ M - |»7s-| = 1} = 


J I[o,b(x,8,T l .-)] (u)N(ds, dx, du, di) 


(0,t] xl d x [0,oo) x r 2 


and 


dt = #{ s '■ \Vs\ ~ \rjs-\ = -1} = / 

J( o 

Note that \r]t\ = bt — dt + |??o| and 6^ are the moments of jumps for q : = bt + dt, so that 


Y I [0Axi,r,Vr-)]( v ) r Vq(dr,d v ). 
' (0,t] X [0,oo) qe xuj, 

XqGrjr — 


Ci = ^ < t}, t > 0. 

fcgN 


For n € N define 


S n ) 


cl = 


/ / [0,6(*, 8 ,» 7 ,_)An] (u)I{\x\ < n}N(ds, dx, du, dj) 


(0,t] xl“x [0,oo)xT2 


+ 


L 


Y / M*i,r,»j P _)An](w)-Z’{|®l < n}j q (dr, dv). 
t] x [0,oo) qeXuJ 
x q €r) r - 


Then 


M, 


(n) 


(x)=c ( t n \x)- [ [ (b(x,s, Va -)An)ds- [ [ (, d(x,s,r, s .)An)ds 

J J x:\x\<n J J x:\x\<n 


is a martingale with respect to By the optional stopping theorem EM q ^ At = 0, hence 


e 1 At 

/ / (b(x, s, 77 s _) A n + d(x, s, ry s _) A n)ds < P{6\ < t} < 1. 

J J x:\x\<n 
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Similarly, 


^771+ 1 At 


OmM 


I x:\x\<n 


(b(x, s, r] s _) A n + d(x, s, ?/ s _) A n)ds < P{O m < t}. 


Consequently, 


t 

/ / (b(x, s,r] s _)/\ n + d(x, s,rj s _) A n)ds 

J J x:\x\<n 
0 

oo 0 ™ +1 

< ^2 E / / (6(x, s, A n + d(x, s, r/ s _) A n)ds 

m=l J J x:\x\<n 
Urn 


< ^2 p {°™ < *} = X] ■ P ‘t c * > 

ra=l ra=l 

Letting n —>• oo, we get by the monotone convergence theorem 


t 


/ / (b(x,s,ri s _) + d(x,s,r] s -))ds 
J J K d 


< -Ect 


Only existing particles may disappear, hence the number of deaths dt satisfies 


d t < b t + 17/o | • 


Thus, 

Sc t < 2 Eb t + S1 7 / 0 1 < 2S|r/t| + S|r/ 0 | < oo. (24) 

Since r/t C rj t a.s., Proposition 3.1 implies (21). 

If follows from the above construction, (22) and (23) that (?/t) is a strong solution to (2). 
Similarly to the proof of Proposition 3.1, we can show by induction on n that equation (2) has a 
unique solution on [0, 9 n \. Namely, each two solutions coincide on [0, 9 n ] a.s. Thus, any solution 
coincides with (?/t) a.s. for all t £ [0, 6 n ]. 

□ 

Remark 3.3. Assume that b and d are time-homogeneous. Let r/o be a non-random initial 
condition, r]o = a, a £ To- The solution of (2) with r/o = a will be denoted as (?/(a, t))t>o- Let 
P a be the push-forward of P under the mapping 

113wa (//( ck , •)) € Dr 0 [0, oo). (25) 

It can be derived from the proof of Proposition 3.2 that, for fixed u £ fl, the unique 
solution is jointly measurable in (t,a). Thus, the family {P a } of probability measures on 
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-Dr 0 [0,oo) is measurable in a. We will often use formulations related to the probability space 
(-Dr 0 [0, oo), &(Dr 0 [0, oo)), P a ); in this case, coordinate mappings will be denoted by rjt, 

r] t (x) = x(t), x £ Z?r 0 [0,oo). 

The processes {rjt)tG[ 0 ,co) and (rj(a, -))re[ 0 ,oo) have the same law (under P a and P , respec¬ 
tively). As one would expect, the family of measures {P a ,a £ To} is a Markov process, or a 
Markov family of probability measures; see Proposition 3.8 below. For a measure /x on To, we 
define 


P,,. = j P a v{da). 

We denote by the expectation under P^. 

Remark 3.4. Let b\,d\ be another pair of birth and death rates, satisfying all conditions 
imposed on b and d. Consider a unique solution (rjt) of (2) with rates b\,d\ instead of b,d, but 
with the same initial condition rjo and all the other underlying structures. If for all ( £ D, 
where D £ ^(Fo(M d )) , &i(-,C) = fr('>C)i di(-,C) = d(-, C), then rj t = r) t for all t < inf{s > 0 : 
r] s (f D} = inf{s > 0 : rj s (f D}. This may be proven in the same way as the theorem above. 

Remark 3.5. Assume that all the conditions of Proposition 3.2 are fulfilled except (5). Then 
we could not claim that (21) holds. However, we would still get a unique solution on [0, oo). We 
are mostly interested in the case of a non-random initial condition, therefore we do not discuss 
the case when (21) is not satisfied. 

Remark 3.6. We solved equation (2) w-wisely. We can deduce from the proof of Proposition 
3.2 that 0 n and z n are measurable functions of rj o and N in the sense that, for example, 6\ = 
F\ (?7o■ A 7 ) a.s. for a measurable function F\ : To x T(R + x R d x M + x T 2 ) —>• M+. As a consequence, 
there is a functional dependence of the solution process and the “input”: the process (r]t)t>o is 
some function of r/o and N. 

The following corollary is a consequence of Proposition 3.2 and Remark 3.6. 

Corollary 3.7. Joint uniqueness in law holds for equation (2) with initial distribution v satis¬ 
fying 



< oo. 


As usually, the Markov property of a solution follows from uniqueness. 


Proposition 3.8. The unique solution o,oo) °f (2) is a Markov process. 
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Proof. Take arbitrary t > to > 0. Consider the equation 


Zt(B)= j I[o,b(x,s,r, s -)](u)I{J I[0,dix,r,nr-)]( v h( d r,dv) =0}N(ds,dx,du,d'y) 


(to,t]xSx[0,oo)xr 2 


+ 


/ !{ / / [0,d(x,r,r ?r _)](«)7(^^) = Q}vt 0 (dx,d'y) +r)t 0 (B), t>t 0 . (26) 

J Jto 


Bx r 2 


The unique solution of (26) is (iJ s )se[to, oo)- As in the proof of Proposition 3.2 we can see that 
(Vs)se [to,oo) i s measurable with respect to the filtration generated by random variables of the form 
N([s, q]xBxU xS), and ?ft 0 ( B ), where B £ ^(R d ), i £ Z, to < s < q, U £ &(R+) H € ^> 4 ^ 2 )- 
Poisson point process have independent increments, hence for any U £ &(Dr 0 [to, oo)) 


F{{vt )te[to,oo) ^ ^ I ^to) P{(Vt)t£[to,oo) ^ ^ I Vt 0 } 


almost surely. □ 

Proof of Proposition 2.7. Let Aq be the image of N under the projection 

(s, x, u, 7 ) i->- (s, x, u). 

The process N\ is a Poisson point process on M + xR d x M + with intensity measure dsdxdu. We 
have 

Vt(B)= J I[o ,b(x,s,rj a -)\(u)I{J IfodfarM^iv^id^dv) = 0}N(ds, dx, du, d'y) 

( 0 ,t]xBx[ 0 ,oo)xr 2 

+ J I{J I[o,d(x,r,r,r-)]i v h( dr , dv ) = 0}f)o{dx,d'y) 

BxT 2 0 

= J I[o,b(x,s,r) s -)]( u )Ni{ds,dx,du) + 770 ( 5 ) 

(0 ,t] xBx[0,oo) 

- / J { x <? G r]r-}I[o,d{x, r , Vr -)\(v)l q {dr,dv). 

gem7( 0)t ] x [ 0oo ) 

By Ito’s formula for F £ ^ 

F(Vt) - F(vo) = J I[o,b(x,s,r, s -)](u){F{vs- Ui) - F(rj s -)}Ni(ds,dx,du) 

(0,4] xB(oj,Rj')x [0,oo) 
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I{x q G T] r _}I [0Ax ^ Vr _ )] (v){F(r] s -\x) - F(i ls _)}'y q (dr,dv). 


+ £ / 

9eXu ^( 0) t] x [ 0iOO ) 


We can write 


]( u ){ F (Vs- Ux) - F(rj s _)}Ni(ds,dx,du) 


(0,t]xB(oj,fijf)x [0,oo) 


+ 


/ 


J b(x,s,r] s -){F(r] s - Ux) - F(rj s _)}dxds 

(0,t]xB(o d ,R F ) 

^[o,6(x, S)f?s _)](w){^(%- Ux) - F(r/ s _)}iVi(ds,dx,du), 


(0,t] xB(o d ,Rp) x [0,oo) 

where N = N — dsdxdu. Since F £ the process 

J I[0,b(x,s, Va -)]( u ){ F iVs- Ux) - -F(r/ s -)}iVi(ds,da:,du) 

(0,t] xB(o,j,fli?)x [0,oo) 

is a martingale by item (vi) of Definition 2.1, see e.g. [IW81, Section 3 of Chapter 2]. Similarly, 

2 / J K € r h-} J [ 0 ,d(x,r,rj r — )]( v ){F(t)s-\x) - F(r] s _)}'i q (dr,dv) 

qeiuj ( 0)f ] x [ 0iOO ) 

can be decomposed into a sum of 

+ ^2 d(x,r,rj r -)]{F(ri s -\x) - F(rj s -)}dr 


(o,t] 


x&rir- 


and a martingale. The desired statement follows. 


□ 


3.1 Coupling 

Here we discuss the coupling of two birth-and-death processes. As a matter of fact, we have 
already used the coupling technique in the proof of Proposition 3.2. 

Consider two equations of the form (2), 


^ (fc)(jB) = / / [0A(-,ne ) )i (v)7(dr ’^ ) = 0} 


( 0 ,t]xSx[ 0 ,oo)xr 2 


xN(ds,dx,du,d'y) + j I{J^I [od ^ xr ^k) )] {v)'y(dr,dv)=0}fi 0 {dx,d'y), k= 1,2. 


(27) 


flxr 2 
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(k) 

Assume that the initial conditions £q ' and the rates bk, dk satisfy the conditions of Propo¬ 
sition 3.2. Let (£f^)te[ 0 ,oo) be the unique strong solutions. 


Proposition 3.9. Assume that almost surely C £q 2 \ and for any two finite configurations 

b\(x, s, rj 1 ) < b2(x, s, r/ 2 ), i£l d ,s>0 ( 28 ) 


1 9 

rj C rj , 


and 


d\(x, s, rj 1 ) > d2(x, s, rj 2 ), x£rj 1 ,s> 0 


Then 


Ct (1) cCt (2) , t€[0,oo). (29) 

We will show by induction that each moment of birth for (^'''j/wo.oo) is a moment of birth for 
(£t 2 ^)te[0,oo) too, and each moment of death for (£f^)te[ 0 ,oo) is a moment of death for (C^)te[0,oo) 
if the dying point is in (£^)te[o,oo)- Moreover, in both cases the birth or the death occurs at 
exactly the same point. Here a moment of birth is a random time at which a new point appears, 
a moment of death is a random time at which a point disappears from the configuration. The 
statement formulated above is in fact equivalent to (29). 

Here we deal only with the base case, the induction step is done in the same way. We have 
nothing to show if n is a moment of a birth of (^ 2 ^)te[o,oo) or a moment of death of (Ct^tefo.oo)- 
Assume that a new point is born for (£t )te[o,oc) at T i> 

The process (£'^)te[o,co) satisfies (2), therefore Ah({x} X {r,} X [ 0 , T it£il-)]) ~ 1 - Since 

A 1 ) _ AH (- A 2 ) _ A 2 ) 

sti- So e so — sti-j 


by (28) 


Ni({x} x {n} x [OAtzpTr,^)]) = 1, 


hence 


(< 2 \ 

The case when T 2 is a moment of death for (Q )te[o.oc) is analyzed analogously. 


□ 
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4 Aggregation model: proofs 


The main idea behind our analysis in this section is to couple the process 0 with another 
birth-and-death process, to which we can apply Lemma 5.1. 

To do so, let us introduce another pair of the birth and death rates, b\,d\, and an initial 
condition £o = Vo LI A, such that bi(x,r/) = di(x,rj) = 0 for x ^ A, di(x,rj) = a'^ for x £ A, 
b\(x,rj) < b(x,rj) for all x,rj, and for some constant c > 0 

J bi(x,rj)dx = c\r] PI A|, ry € To- 

A 

There exists a function b\ satisfying these assumptions. 

Functions b\, d\ satisfy conditions of Theorem 2.6. Furthermore, the conditions of Proposi¬ 
tion 3.9 are fulfilled here: for yy^jy 2 £ To, ry 1 C ?y 2 we have 

6i(x, ry 1 ) < b(x, jy 1 ) < b(x,r/ 2 ) 


as well as 

di(x,rj 1 ) > d(x,ri l ) > d(x,r] 2 ). 

Denote by (£t)t> 0 the unique solution of (2) with the birth and death rates bi,d± and initial 
condition £o- By Proposition 3.9, C ryj hold a.s. for all t > 0. 

In this section we will work on the canonical probability space 

(Dr 0 [0,oo) x D r0 [0,oo),^(Dr 0 [0,oo) x D To [0,oo)),P Q ), 

where P a is the push-forward of the measure P under 

P 3 w (ry(a, •), (£(a, •)) € D r 0 [0, oo) x D Vo [0, oo). 

Consider the embedded Markov chain of the process (£t)t>o, hfc := £,r k , where are the 
moments of jumps of (£(). It turns out that the process u = {u^jfceN) where Uk ■= |Yfc|, is a 
Markov chain too. Indeed, the equality 

P ai {\Yi\ = k} = P a2 {\Yi\ = k}, ken,ae T 0 . 
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holds when \a\ fl A| = |«2 H A|, since both sides are equal to 


_C_ 

c+a — l Q! i nA l 

a -|g 1 nA| 

c+a _ l Q: i nA l 

0 


if k = |«i n A| + 1, 
if k = \a\ PI A| — 1, 
in other cases. 


Therefore, Lemma 5.1 is applicable here, with /(•) = | • |. 

Proof of Proposition 2.8. Having in mind the inclusion ^ C tjt (P a -a.s.), we will prove 
this lemma for (£t). 

The transition probabilities for the Markov chain {iik}k&z + are given by 


Pi,j = Pa{uk = j | \uk-l \ =l} = 

for i £ N, j £ Z_|_, and po.j = /{ J= o}, see Remark 5.2 and (6). 

Since the zero is a trap and it is accessible from all other states, there are no recurrent states 
except zero, and the process u has only two possible types of behavior on infinity: 


c+a * 
q~ i 

c+a~ z 

0 


if j = i + 1, 
if j = i - 1, 
in other cases, 


(30) 


P q { 3/ £ N s.t. ui = 0 or lim u m = oo} = 1. 

m—Kx) 


We will now use properties of countable state space Markov chains, see, e.g., [Chu67, § 
12, chapter 1], Chung considers there Markov chain with a reflecting barrier at 0, but we 

m 

may still apply those results, adapting them correspondingly. Denote Q m = W • Then 


k =l 


the probability P a {3k € N s.t. u & = 0} equals to 1 if and only if ^2 Qj = oo, whichever initial 

3 = 1 

OO 

condition a, |aDA| > 0, we have. Moreover, if ^ Qj < oo and P Q {uo = q} = 1 (or, equivalently, 

3= 1 


|a fl A| = q), then p q := P Q {3A: £ N s.t. = 0} = 


E Qj 


■ j(j+l) , 
Qj = c J a 2 ; and 


!+ E Qj 

3 = 1 


From (30) we see that in our case 


Pq = 


~ . 3U+1 ) 

c J a 2 
1 =<?_ 

, “ • j(j+D 

1+^0 c _ ^a 2 
l=i 


< 


E c j a _ 
j=q 


(31) 


i + E c ~^ a 2 
1=1 


Now, for arbitrary C > 1 choose q £ N for which c x a 2 < C 1 . For j > q we have 
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^ 3 n 2 


a 2 < c J a 


-7 ~ 2R 

J n 2 


(c l a 2 q y < C 3 . and 


E 


c 33 a 


j=q 


<Y,c-’ = 


3=1 


c q 

i -c- 1 ’ 


so that the statement of the lemma for (£t)t >0 follows from (31). □ 

Note that the number of particles of the process will go to infinity with probability 1 even 
though the probability of extinction is positive, unless £>(•, 0) = 0 almost everywhere with respect 
to the Lebesgue measure. However, if £>(•, 0) = 0, then 


P ( {\£t\ = 0 for large t}U{|£t |—> oo ,t —> oo} ) = 1 


P ( {|£t| = 0 for large t}n{|£i|—»oo,t—>■ oo} I = 0. 


and 


The following equality is also taken from [Chu67, § 12, chapter 1]; for q > s and all /3 with 
1/3 n A| = q, 


E Qj( s ) 

Pp{3k € N : \u k \ = s} =-=^-, 

1 + E Qj( s ) 

j=s+l 


m , 

where g m (s) = FT = c -(m-s)- 1 (m-s)(m+s+i) in our case 

V , , , Pk,k+1 

k=s +1 


J2 c -0'- s ) a _ 3^-s)(i+s+i) 

P^eN:|ii fc | = s} =---- :=Cg,*<l- (32) 

1+ ^ c -0-s) o -20-s)0+s+i) 

j=s+l 


Note that 


Cg+ 1,1 -t 0, q -)■ oo (33) 

Proof of Proposition 2.9. Let (Xfc)fc G ^ be the embedded chain of (r)t)t> 0 - First we will 
show that for all m G N and a £ To, 

P a {\X k n A| = m infinitely often } = 0. (34) 

Let (3 G Fq, |/3 (~l A| = m, m € N (the case of m = 0 is similar, and we do not write it down). 
Denote k = min{/c € N : X k nA^ Xq n A}. Since C r) t holds - a.s., 


23 






Pp{\Xk n A| > m, Vfc > k } > Pp {(Yfc n A| > m,Mk > l} 
= Pgjufc > rri, Vk > l}. 


( 35 ) 


By (32), the probability Pp{iik > m,Vfc > 1} is positive and does not depend on ( 3 , |/ 3 D A| = 
m: 

®ra • — Pf}\y>k P m,Vk ^ 1} P Pm,m- |-l(l Cm-|-l,m) O' (36) 

Define fcf 1 , i £ N, subsequently by = min {A: > k'j 1 : \X], n A| = m and 3k < k : 
\X^ n A| / m}, k™ = 0. Note that for all /? 


P^Ebio : \X n n A| = m for all n > no| = 0. 


By the strong Markov property, 


P a < \Xk D A| = m inhnitely often > < P a < k 1 ] 1 < oo, Vj G N 


.0 Pa{k'j+ i < oc \ kj 1 < oo} = 0, 

j= 1 


by (35) and (36). Indeed, if P a {k^ < oo} > 0, then 


P a {kj\_ i < oo | kj 1 < oo} = 


P&I{k r ! rl <ioo}PX/ c m{k-i ^ Oo} 


E n L 


Q 4{fc’7 1 <oo} 


< 


E a I{k™< oo}(! - p x k rn{ \X k n A| > m,Vk> k }) 


E n I, 


OL J -{k 7 7 l < oo} 


^^-< 00 }(! - Px k m{u k > m,\/k > 1}) 

<----= 1 - s m < 1 . 

-^ / a^{kj l < 00 } 


Having proved (34), we observe that 


{\Vt n A| — > 00 } U {dt' : Vt > t ', \rjt fi A| = 0 } 

OO 

[^J {\Xk D A| =m infinitely often} 

m= 1 

Note that if for some element of probability space u € D the process (r]t)t >0 is stuck in a trap 
7 , 7 fl A = 0 , then 00 belongs to the set on the left-hand side of (37) and does not belong to the 
set {\Xk fl A| = m inhnitely often}, m G N. 
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The statement of the proposition follows from (34) and (37). □ 

Proof of Proposition 2.10. Define rjt ■= rjt fl A and let X & = rjj k , where ft is the ordered 
sequence of jumps of (jjt)t> 0 - Of course, the process {rjt}t>o is not Markov in general, and neither 
is {**}* e ^. However, for all a G To(M d ) the inequality 


P a{\Xi\ - |* 0 | = !} > P|anA|,|anA|+l 

holds, because for every £ G To, ( PI A = m, the integral of the birth rate &(-, over A is larger 

than cm, and the cumulative death rate in A, d(x,(), is less than ma~ m . 

zeCnA 

The probability of the event that absolutely no death occurs is positive, even when the initial 
configuration contains only one point inside A: 


P a \ I hi I — | hi— | > 0 for all t > 0 > = P a < |*fc+i| — |**| = 1 for all k G N 


ll P a {\X k+1 \ - |* fc | = 11 |**| - |**_i| = 1, |*! | - |*01 = l} 


fceN 


all mihi-iv„i = i} 


fce n 


Cer°(R d ) 
|CnA|=|anA|+fc 


OO 

an 


Pi,i +1 — 


n 


H“l 


c + a~ l 

i=\a\ i=\a\ 


n (i - —— o > o, 

11 v c + a~ lJ 


because the series c +a- i conver g es - In particular, }} Pi^+i —* 1 as m goes to oo. Also, 


i=\a\ 


p a n {\r]t\ - |hi— I > 0 for all t> 0} —>• 1, \a n fl A| -S- oo. 


(38) 


It is clear only an a.s. finite number of deaths inside A occurs on \3t' : Vt > t' , |rjt D A| = 0}. 
By Proposition 2.9, it remains to show that only an a.s. finite number of deaths inside A 
occurs on {| r] t fl A| —>• oo} = {\rj t \ —> oo}. Let us introduce the stopping times a n = inf{s G R : 
\rj s \ > n}, which are hnite on {\rjt\ —> oo}- Only a finite number of events (births and deaths) 
occur until arbitrary hnite time Pp- a.s. for all /? G To, hence for n G N 

Pa(^{\vt\ ~ |hi— I A 0 for all but finitely many t > 0} D {\rjt\ ~> oo}j 
> P a ({\rj t \ - I hi— I > 0 for all t > <r n } n {\rj t \ -s- oo}) 

= Palm-^ojP^m ~ lht-1 > 0 for all t > 0}. 
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From \ri an \ > n we have by (38) 

P Va n {\Vt\ - \rjt-\ > 0 for all t > 0} -7- 1, n -7 oo. 


Therefore, 


P 0 (y{Vh\ — \vt-\ > 0 for all but finitely many t > 0} D {\r] t \ —> oo}j = P a {\r]t\ — > oo}. 


□ 

Proposition 2.10 is also applicable to (£)t> o, since b\,d\ satisfy all the conditions imposed 
on b, d. 

Proof of Theorem 2.11. First we prove the Lemma for (£)t> o'. we prove that for P a -almost 
all oj E F\ := { lim FI A| = oo}, 

t—> OO 

lim inf Kill—i > 0 . ( 39 ) 

t —^oo e ct 

There is no loss in generality in assuming uq = |an A| > 0. Let 0 = tq < t± < T 2 < ... be the 
moments of jumps of (£t)t>o, so that £ rj , = Y k . We recall that the random variables u n = \Y n \ 
constitute a Markov chain. Denote = cn + na~ n . Then 

/ bi(x,Y k )dx + ^2 d 1 (x,Y k ) = c\Y k \ + \Y k \a~^ = ip(u k ). 
a seyjb 

By Theorem 12.17 in [Kal02] the random variables ^ k := V’( w fc)( r fc +1 — Tfc), fc E Z + are 
independent and exponentially distributed with parameter 1. Furthermore, the sequence { 7 ^} 
is independent of Y. In particular, it is independent of {u k } k£ z+- 

From Proposition 2.10 we know that only a finite number of deaths inside A occur a.s. on 
F\ . Particularly, there exists a positive finite random variable m such that the inequalities 

uq + n> u n > uq + n — m(w), n E N (40) 


hold with probability 1. 

We can write 

7i—l n— 1 n—1 

r„ = Efe+i * 7 = E ^ E 

k =1 k =1 rv k =1 u 

00 

Due to Kolmogorov’s two-series theorem, the series ^ U Q+ck divergent (we recall that 

/c=l 

E^ k = = 1). Hence r n —> 00 a.s. 
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We will show below that 


(41) 


cr n < In n + cry, n G N, 

where 7 is some finite random variable. Using (41), we obtain 


0 ct 


- ( m+ l) e C7 jt - 0 } = Pa {\^n, \ > (rrj G N ) 


(m + l)e c 7 ’ 


= P a {u n >- -e CTn+1 c y n G N} = P a { In (u n ) + ln(m + 1) > cr n+1 - Py,n € N} = 1. 

k m + 1 ’ K ’ 

Therefore, (39) holds. 

Inequality (41) follows from the convergence of the series 


£ 


Ik 1 


ck J 

To establish the convergence of (42), we note that 


(42) 


£ 


Ik Ik 


—[ \^{Uk) cu k 


converges by Kolmogorov’s theorem: 


(43) 


u k 


sr' ( ik _ 7 fc_\ — w Uka Uk < _L 1 

“ \ip(uk) cu k ) ^^ k cu k ip(u k ) _ C 2 Uk 


^ rn ^ 00 


1 m — it 1 -t 00 

1 \ - a k 1 \ - a 


-3 


and 


1 1 y—-r l Li l y^ a 

k=1 /c=m+l k=1 j'=l 


vfl- — 

“ V cfc 


(44) 


converges by (40). 

The convergence of the series in (42) follows from the fact that (43) and (44) converge. 

We have thus proved (39). To establish the statement of the theorem, note that a n = inf{f > 
0 : |? 7 t| > n} is finite on F and 

{ liminf ^ = 0, \rjt \ —> 00 } C { liminf ^ = 0}. 


It follows from what we have already proved that 


16 


Pp{ liminf = 0 } = Pp{{£,t)t >0 extincts}, /3 G T 0 . 
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Therefore, by Proposition 2.8 and the strong Markov property (for the latter see e.g. [Kal02, 
Theorem 12.14]) 

P a { liminf ' A = 0, |? 7 f | 00 } = P a P,- n { liminf = 0, |? 7 t | -> 00 } 

<P«P,»„{liminf^!=0}<C-", 

where C is the constant that appeared in Proposition 2.8. Since n is arbitrary, 

P Q { liminf ^ = 0, \rj t \ ->■ 00 } = 0. 


□ 

Proof of Corollary 2.12. Let us fix a configuration a, aflA/0 , We saw in the proof of 
Theorem 2.11 that for almost all cj £ F = {cj : liminf > 0} we have 

t —>00 e 


p a m > 


1 


(m + l)e cX < 


-e , t > 0} = 1. 


Let F k be the set {u : h ^ > £}. Then ]J F k = F, and, since P a (F ) > 0, 

v ’ ke n 


Pa{F k ) > 0 


for some k € N. Hence 


E a \rj t n A| > E a \r] t D A| I Fk > ^e ct P a (F k ). □ 


5 Appendix 

5.1 Markovian functions of a Markov chain 

Let (S,&(S)) be a Polish (state) space. Consider a (time-homogeneous) Markov chain on 
(S,&(S)) as a family of probability measures on S°°. Namely, on the measurable space 
(fl,J^) = (S°°, consider a family of probability measures {-P s } sG s such that for the 

coordinate mappings 


X n :Q^S, 

X n (s±, S2,...) — S n 

the process X = {X n } ne z is a Markov chain, and for all s € S 
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Ps{X o = s} = 1 


Ps{ Ara+mj S Aj, j — 1, k\ | t^n} — Px n {X m j S Aj,j — 1, —, /ci}. 

Here Aj £ £§{S), nij £ N, k± £ N, = a{X\,..., X n }. The space S is separable, hence 
there exists a transition probability kernel Q : S x BS{S) —>• [0,1] such that 

Q{s,A) = PsiXx £ A}, s £ S, A £ @(S). 

Consider a transformation of the chain X, Y n = f(X n ), where f : S —> Z + is a Borel- 
measurable function, with convention ^(Z+) = 2 Z +. In this section we will give sufficient 
conditions for Y = {Y n } ne z, + to be a Markov chain. A similar question was discussed in [BR58]. 

Lemma 5.1. Assume that for any bounded Borel function h : S —» S 

E s h{X 1 ) = E q h(X 1 ) whenever f{s) = f(q), (45) 

Then Y is a Markov chain. 

Remark. Condition (45) is the equality of distributions of X\ under two different measures, 
P s and P q . 

Proof. For the natural filtrations of the processes X and Y we have an inclusion 

n£ N, (46) 

since Y is a function of X. For k £ N and bounded Borel functions hj : Z + — > M, j = 1, 2,..., k 
(any function on Z + is a Borel function), 


E„ 


II hj(Y n+j ) I 


3 = 1 


Ex n n hjifiXj)) 


3 =1 


/ Q(xo,dxi)hi(f(xi)) / Q(xi,dx 2 )h 2 (f{x 2 ))... / Q(x n -i,dx n )h n (f(x n )) 
Js Js Js 


(47) 


Xq — X n 


To transform the last integral, we introduce a new kernel: for y £ f(S) chose x £ S with 
f(x) = y, and then for B C Z + define 


Q(y,B) = Q(x,f~ 1 (B)). 


(48) 


The expression on the right-hand side does not depend on the choice of x because of (45). To 
make the kernel Q defined on Z + x <S$'(Z+), we set 
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Q{y,B) = /{ogb}> v i fi s )- 


Then, setting z n = f(x n ), it follows from the change of variables formula for the Lebesgue 
integral that 


I Q(x n — i, dxn^hn (/(Xfj,)) — / Q(f (x n —i), dz n ')h n (z n '). 

Js JZ+ 

Likewise, setting z n -\ = f(x n - 1 ), we get 

/ Q(x n — 2 idx n —\)h n ( y f{x rL —\)) / Q(x n — i,d.'r n )/in,(/(x n )) — 

Js Js 

/ Q{x n — 2 idx n —\)h n ( K f{x n —\)^j / Q(yf(yX n —\)i dz n )h n ( y z n ') — 

JS Jz + 

/ Q(f(%n— 2 ); dz n — i)/i n (z n _ 1 ) / Q(^z n —i, dz n ')h n (^z n '). 

J Z+ J Z_|_ 

Further proceeding, we obtain 

/ Q{x 0 ,dxi)hi(f{xi)) / Q{xi,dx 2 )h 2 {f{x 2 ))... / Q{x n - 1 ,dx n )h n (f(x n )) 

Js Js Js 


Q(zo,dzi)hi(zi) 


Q{zi,dz 2 )h 2 {z2)— / Q(2„-i,d2 re )/i n (^ re ), 


where z 0 = f(x 0 ). 
Thus, 




i=l 


/ Q(f(X 0 ),dzi)hi(z 1 ) / Q(zi,dz 2 )h 2 (z 2 )... / Q(z n - 1 ,dz n )h n (z n ). 

J Z+ J Z+ J Z-j- 

This equality and (46) imply that L" is a Markov chain. 

Remark 5.2. The kernel Q and the chain f(X n ) are related: for all s € S, n, m € N and 
M C N, 


P,{/(X n+1 ) € M | /(X n ) = m} = Q(m,M) 

whenever P s {/(X„ + i) = m} > 0. Informally, we can say that Q is the transition probability 
kernel for the chain {/(X n )} ng g + . 

Remark 5.3. Clearly, this result holds for a Markov chain which is not necessarily defined on 
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a canonical state space, because the property of a process to be a Markov chain depends on its 
distribution only. 

5.2 Strong Markov property of the driving process 

Let N be compatible with a right-continuous complete filtration {^t}t> o> and r be a finite a.s. 
{^i}t>o-stopping time. For 7 G r 2 , 7 = ^ 8^ Ui) , let (9 r 7 = ^ <5 (s ._ T) „.). Also, for S € ^(r 2 ) 

i i:Si>T 

we define the shift 

0 T H = {7 6 L 2 | 6» r y € H} 

Introduce another point process A" on x R d x x T 2 , 

A([0, s] x U x S) = N((t, t + -s] x U x 0 T S), s > 0, 17 G ^(M d x R+), E € #(r 2 ). 

Proposition 5.4. T/ie process N is a Poisson point process with intensity measure ds x dx x 
du x 7r, independent of 3P r . 

Proof. To prove the proposition, it is enough to show that 

(i) for any b > a > 0, open bounded U C M d X M + and open T C r 2 , IV((a, b) X U x E) is a 

Poisson random variable with mean (b — a) X l(U) x 7r(S), where l is the Lebesgue measure on 
the corresponding space, and 

(ii) for any bk > a k >0, k = 1 open bounded Uk C K d and open S& C T 2 such that 

((a*, bi) xUiX 5 i) n (( aj,bj ) x Uj x S f) = 0, i ± j, the collection {N((a k , b k ) xU k x ’E k )}k=i,m 
is a finite sequence of independent random variables, independent of & T . 

Let r n be the sequence of {^i}t>o-stopping times, r n = Jr on {r € (^r, Jr]}, k G N. Then 
T n f t and T n — t < Jr. The stopping times r n take only countably many values. Therefore the 
process N satisfies the strong Markov property for r n : the processes N n . defined by 

A n ([0, s\ x U x E) := N((r n , r n + s] x U x 0 Tn 5), 

are Poisson point processes, independent of J^ Tn . 

To prove (i), note that N n ((a , b) xU x E) — > N((a, b) xU x 5) a.s. and all random variables 
N n ((a, b) x U x E) have the same distribution, therefore A((o, b) x U x E) is a Poisson random 
variable with mean (6 — a)l(17)7r(S). The random variables N n ((a,b ) X U x E) are independent 

of PP T . : hence N((a,b) x 17 X E) is independent of d^, too. Similarly, (ii) follows. □ 
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